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It is well known that the late-time behaviour of gravitational collapse is 
dominated by an inverse power-law decaying tail. We calculate higher-order 



tional collapse. The dominant "contamination" is shown to die off at late 



CN ' corrections to this power-law behaviour in a spherically symmetric gravita- 
> 

o 

' times as M^t '^ln(t/M). This decay rate is much s/ower than has been consid- 

o . 

' ered so far. It implies, for instance, that an 'exact' (numerical) determination 

^ ■ of the power index to within ~ 1% requires extremely long integration times 
I ' 

^J^' of order lO^M. We show that the leading order fingerprint of the black-hole 



2+-4 



electric charge is of order Q t 

I. INTRODUCTION 

Waves propagating in a curved spacetime develop "tails". In particular, it is well es- 
tablished that the dominant late-time behaviour of massless fields propagating in black-hole 
spacetimes is a power-law tail. Price [|I| was the first to analyze the mechanism by which 
the spacetime outside a (nearly spherical) collapsing star divests itself of all radiative mul- 
tipole moments, and leaves behind a Schwarzschild black hole; it was demonstrated that 
all radiative perturbations decay asymptotically as an inverse power of time. Physically, 
these inverse power-law tails are associated with the backscattering of waves off the effective 
curvature potential at asymptotically far regions pllll. 



The analysis of Price was extended by other authors. Bicak generalized the analysis 
and studied the dynamics of a scalar field in a charged Reissner-Nordstrom spacetime. He 
also found an asymptotic inverse power-law decay of the field, with the same power indices 
as in the Schwarzschild spacetime (with the exception of the extremal Reissner-Nordstrom 
black hole [^). In a brilliant work, Leaver |^ demonstrated that the late-time tail can be 
associated with the existence of a branch cut in the Green's function for the wave propagation 
problem. 

Gundlach, Price, and Pullin |^ showed that these inverse power-law tails also charac- 
terize the late-time evolution of radiative fields at future null infinity, and at the black-hole 
outer horizon. Furthermore, they showed that power-law tails are a genuine feature of grav- 
itational collapse - the existence of these tails was demonstrated in full non-linear numerical 
simulations of the spherically symmetric collapse of a self-gravitating scalar field (this 
was later reproduced in [0]). 

Our current understanding of the late-time tail is, however, somewhat unsatisfactory. 
The {leading order) power-law tails in black- hole spacetimes are well established ]I|-pi?|, 
but the resultant formulae are only truly useful at very late times. In a typical evolution 
scenario there is a considerable time window in which the signal is no longer dominated by 
the quasinormal modes [^|Tli^|6|,|T2|,|8|,|9| , but the leading order power-law tail has not yet 
taken over 0. 

The purpose of this paper is to derive analytic expressions for the higher-order corrections 
which "contaminate" the well-known power-law tail in a spherically symmetric gravitational 
collapse. The determination of these higher-order terms is important from several points 
of view: The analyses of Bicak and Gundlach et. al. established the fact that the 
leading-order power-law tail is universal in the sense that it is independent of the black-hole 
electric charge (i.e., the power index in a c/ianjfec? Reissner-Nordstrom spacetime was shown 
to be identical with the one found by Price for the neutral Schwarzschild black hole). 
This observation begs the question: what fingerprints (if any) does the black-hole electric 
charge leave on the field's decay ? 
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Moreover, the calculation of higher-order corrections to the leading order power-law tail 
is also of practical importance; this is especially crucial for the determination of the power 
index from numerical simulations. The dominant inverse power-law tail is "contaminated" hj 
higher-order terms, whose effect become larger as the aveliable time of integration decreases. 
The precise power index is expected only at infinitely-late time. Thus, in practice, the limited 
time of integration introduces an inherent error in the determination of the power index. The 
only systematic approach to quantify the errors which are introduced by the finite integration 
time is to study higher- order corrections. If one computes the contaminated part of the 
late-time tail, then the ratio of the corrections to the leading order term is a systematic, 
quantitative, indication of the error caused by the finite-time numerical calculation. 

These questions and several others are addressed in the present paper. The plan of 
the paper is as follows. In Sec. we give a short description of the physical system and 
formulate the evolution equation considered. In Sec. |T| we give an analytical description of 
the late-time evolution of scalar fields in black-hole spacetimes. In Sec. ^ we confirm our 
analytical results by numerical simulations. We conclude in Sec. [V] with a brief summary of 
our results and their implications. 



II. DESCRIPTION OF THE SYSTEM 

We consider the evolution of a spherically symmetric massless scalar field in a spherically 
symmetric charged background (a collapsing star or a fixed black hole). The external grav- 
itational field of a spherically symmetric charged object of mass M and charge Q is given 
by the Reissner-Nordstrom metric 

ds^ = -[l + \\dt^+il + ^ dr^ + r^dn^ . 1 

Using the tortoise radial coordinate y, which is defined by dy = dr/{l — 2M/r + (5^/r^), the 
line element becomes 

ds^=(^-'^-^ + ^-de + dy^) + T^dn^ , (2) 
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where r = r{y). 

The wave equation 4'-ab9°'^4' = for the scalar field in the black-hole background is 

i^,tt - ^,yy + \/V = , (3) 

where 

^ = WO=(l- — + -^j • (4) 

In terms of the tortoise coordinate y and ioi y ^ M the curvature potential Eq. (1) reads 

v=v^,^ = ^J^^llM^.^Ml±^^o(^]. (5) 

yi yi yA ^ y5 J 

III. THE LATE-TIME BEHAVIOUR 

The general solution to the wave-equation can be written as a series depending on 
two arbitrary functions F and G [|I|] 

oo 

+ j:Bk{y) [G^-'-'\u) + {-lfF^-'-'\v)] . (6) 

fc=0 

Here u = t — y is a retarded time coordinate and v = t + y is an advanced time coordinate. 
For any function H, H'^^^ is the fcth derivative of negative-order derivatives are to be 
interpreted as integrals. The first two terms in Eq. (|^) represent the zeroth-order solution 
(with V = {)). 

The star begins to collapse at a retarded time u = uq. The world line of the stellar surface 
is asymptotic to an ingoing null line v = Vq, while the variation of the field on the stellar 
surface is asymptotically infinitely redshifted This effect is caused by the time dilation 
between static frames and infalling frames. A static external observer sees all processes on 
the stellar surface become "frozen" as the star approaches the horizon. Thus, he sees all 
physical quantities approach a constant. We therefore make the explicit assumption that 
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after some retarded time u = ui, du4> = on f = iiq. This assumption has been proven to 
be very successful p|,|5|,pr3[] . 



We begin with the first stage of the evolution, i.e., the scattering of the field in the region 
Uo < u < Ui. The first two terms in Eq. (|^) represent the primary waves in the wave front, 
while the sum represents backscattered waves. The interpretation of these integral terms as 
backscatter comes from the fact that they depend on data spread out over a section of the 
past light cone, while outgoing waves depend only on data at a fixed u [Q. 

After the passage of the primary waves there is no outgoing radiation for u > ui, aside 
from backscattered waves. This means that G{ui) = 0. Hence, for a large y a.t u = ui, the 
dominant term in Eq. (H) is ^{u = Ui,y) = BQ(y)G^^^\ui). 

The functions Bk{r) satisfy the recursion relation 2B'^ — B'1_^ + V Bk_i = for A; > 1, 
where B' = dB/dy^ and B'q = —V{y)/2. Thus, one finds 

^(M = ni,y) = [My-V2 + 2MV'lny-QVV3]G'(-i)(«i)[l + 0(M/y)] . (7) 

This is the dominant backscatter of the primary waves. 

With this specification of characteristic data on u = ui, we shall next consider the 
asymptotic evolution of the field. We confine our attention to the region u > ui, y ^ M, \Q\. 
In this region the spacetime is approximated as flat [|l],H . (The validity of this approximation 
is ultimately justified by numerical simulations). Thus, the solution for ip can be written as 

^ = ^(o)(«) + /(o)(^) . (8) 
Comparing Eq. (|) with the initial data on u = Ui Eq. (|^), one finds 

f{v) = FqV-^ + Fiv'^ \nv + F2V-^ , (9) 

where 

Fo = 2MG^'^\ui) , Fi = 16M^G^-^\u^) , Fa = -8Q^G^-^\u^)/3 . (10) 

oo 

For late times t y we can expand g{u) = {—l)"g^"''(t)y"'/n\ and similarly for /(f). 

n=0 

Using these expansions we can rewrite Eq. (|^) as 



n=0 



(11) 



where the coefficients Kq are those given in 

Using the boundary conditions for small r (regularity as y ^ —oo, at the horizon of 
a black hole, or at r = for a stellar model), one finds that at late times the terms 
h{t) = f{t) + g{t) and f^^\t) must be of the same order (see for additional details). 
Thus, we conclude that 



f{t) ^ For^ + FirHnt + F2t~^ , g{t) ^ - Fot"^ + Fi^^lnt 



(12) 



and 



h{t)^f{t) + g{t) = 0{MH-') 



(13) 



We therefore find that the late-time behaviour of the field for t ^ y ^ M, \Q\is 



= -2Kly \2Fot-^ + SFit"^ ln(t/M) + SFat"^! fl + 0{M/t) 



(14) 



This is the late-time behaviour of the field at a fixed radius. 



IV. NUMERICAL RESULTS 



It is straightforward to integrate Eq. (^) using the method described in |5|. We have 
used, however, a modified version of the numerical code used in |T^, which is essential to 
achieve the extremely high accuracy needed for the computation (see for additional 
details). 

The late-time evolution of the scalar field is independent of the form of the initial data 
used. The results presented here are for a Gaussian pulse on u = 

'4,{u = ^,v) = Ae^^[-[{v-vo)laf] , (15) 
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with a center at vq = 10 and a width cr = 2. The black-hole mass is set equal to M = 0.5; 
this corresponds to the freedom to rescale the coordinates by an overall length scale. 

The temporal evolution of the field at a fixed radius y = 50 is shown in the top panel of 
Fig. |I[ The dominant power-law fall off is manifest at asymptotic late times. In order to 
study the contamination effect of higher-order terms [see Eq. (Ill)], we use the notion of a 



local power index 'j, defined by 7 = —tipt/ip [0]. Taking cognizance of Eq. (Q) we find 

M 

7 = 3 + 12— ln(t/M) + 0{M/t) . (16) 

The approach of the local power index to its well-known asymptotic value 7asj, — > 3 is depicts 
in the bottom panel of Fig. |l|. The plot shows that 7^3 from above, with a qualitative 
agreement with Eq. (p^Gf). 



In order to establish quantitatively the physical picture presented in Sec. |T|, we define 
the quantity 5 = (7 — 7asy)(t/M)/ ln(t/M). Figure 0. depicts 5 as a function of t at three 
surfaces of constant radius y = 5, 10, and 50 (from bottom to top). The numerical result 
Sasy — 12 (independently of the value of y) is in excellent agreement with the analytically 
predicted behaviour 6 ^ 12 [see Eq. (0)]. Thus, Fig. 0. establishes the existence of the 
contamination term of order M^t~^ In (t/M). It should be noted that the value of 'jasy used 
in this figure is 2.998 rather than the theoretical value 3. This slight deviation from the 
theoretical value is expected due to the corrections of order 0{M/t) in the expression for 7 
Eq. m. 



V. SUMMARY AND PHYSICAL IMPLICATIONS 



The purpose of this paper was to derive analytic expressions for the higher-order correc- 
tions which "contaminate" the well-known power-law tail in a spherically symmetric grav- 
itational collapse. We have shown, both analytically and numerically, that the dominant 
correction dies off at late times as M^t~^ ln(t/M). This late-time decay of the contamination 
is much slower than has been considered so far ||^ (see the discussion in Appendix B) . 
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Aside from being theoretically important, the result Eq. (JT^ is also of practical impor- 
tance. It follows that an 'exact' (numerical) determination of the power index demands 
extremely long integration times. The most accurate method for determining the power in- 
dex experimentally applies to the concept of the local power index. In this way one discard 
the relatively large contamination which characterizes the early stages of the evolution. Still, 
it follows from Eq. (^) that a determination of the power index to within ~ 1% requires 
an integration time of order t ^ 3000M . 

The dominant power-law tail is known to be universal in the sense that it is independent 
of the black-hole parameters 0J^. We have shown, however, that this universality is removed 
once we consider higher-order corrections terms - the leading order fingerprint of the black- 
hole electric charge behaves as Q'^t~'^. 
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APPENDIX A: LATE-TIME BEHAVIOUR AT NULL INFINITY AND AT THE 

HORIZON 

It follows from Eqs. (||), (P), and (|12D that the asymptotic behaviour of the field at 
future null infinity scrz+ (i.e., at v ^ u) is 

V'(t^ > ^ 5''^°^") - - (-^oM~^ + -^iW^lnw) . (Al) 

We finally consider the behaviour of the field at the black- hole outer horizon r+. As 
y — cxD the curvature potential Eq. (^ is exponentially small, and the general solution 
to Eq. (|^) can be written as ip = a{u) + 7(f). On t> = t>o we take ducf) = (for u 00). 
Thus, a{u) must be a constant, and with no loss of generality we can choose it to be zero. 
We next expand 7(f ) for t |y| as 
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OO -| 

^ = ^(^) = ^^W(t)^«. (A2) 



In order to match the y <^ —M solution Eq. ( [A^ ) with the y ^ M solution Eq. (|14|) 



we make the ansatz ~ ipstatiu) + 3Fit~^ ln(t/M) + 3F2t~'^] for the solution in the 

region y <^ —M and t ^ \y\. In other words, we assume that the solution in the y <C —M 
region has the same late time t dependence as the y ^ M solution. This assumption has been 
proven to be very successful for the leading order behaviour of both neutral and charged 
g fields. Using this assumption, one finds -f{t) = Tq [2Fot-^ + 3Fit-^\n{t/M) + SFzt"^], 



where Fq is a constant. Thus, the asymptotic behaviour of the field at the black- hole horizon 
is 

^(m ^oo,v) = To \2Fov-^ + SF^v'^ ln{v/M) + 3F2V~^] . (A3) 



APPENDIX B: A NOTE ON ANDERSSON'S ANALYSIS 

The first attempt to calculate higher-order corrections to the dominant power-law tail 
was made by Andersson (in the context of the Schwarzschild spacetime). This analysis 
was based on an approximated curvature potential in the region far away from the black hole 
[see Eq. (26) of P]]. While this approximated potential simplifies the analysis, it actually 
misses the genuine leading order corrections terms. The leading order correction in [P] was 
found to be of order Mt~^. This (artificial) result is caused by the fact that terms of order 
M/ (and smaller) were totally neglected by the approximated curvature potential used in 

i- 

Moreover, the approximated approach presented in 0, z/ extended along the same lines 
to a Reissner-Nordstrom spacetime would imply that the influence of the black-hole electric 
charge on the late-time tail vanishes identically (i.e., it would vanish to any order in t~^). 
This result is again a direct consequence of the fact that the approximated approach of 
does not take into account curvature terms of order Q'^/r^ (and M'^/r'^) which appear in the 
(exact) curvature potential. 
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FIG. 1. Temporal evolution of the scalar field, evaluated at y = 50 in a Schwarzschild space- 
time with M = 0.5. The initial data is a Gaussian distribution with vq = 10 and a = 2. The 
asymptotic power-law fall off is manifest at late times (top panel). The bottom panel depicts the 
evolution of the local power index 7 = —ttpt/'^p- The power index approaches the well-known 
asymptotic value ^yasy = 3 at late times. 
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FIG. 2. The time evolution of the quantity ^ = (7 — jasy){t/M)/ln{t/M), evaluated at 
y = 5, 10, and 50 (from bottom to top). The asymptotic value Sasy — 12 is in excellent agreement 
with the analytically predicted behaviour (5 — *■ 12 at late times. This result establishes the existence 
of the contamination term of order M^t~^ ln(i/M ). 
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